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Abstract. We investigate the relationship between the C"*0 column density and the visual extinction in Chamaeleon I and in 
a part of the Chamaeleon III molecular cloud. The C"*0 column densities, A'(C'*0), are calculated from 7=1-0 rotational 
line data observed with the SEST telescope. The visual extinctions, Ay, are derived using JHK photometry from the 2MASS 
survey and the NICER color excess technique. In contrast with the previous results of Hayakawa et al. ( .2001,) , we find that 
the average A?(C'^0)/Av ratios are similar in Cha I and Cha III, and lie close to values derived for other clouds, i.e. A'(C'^O) 
» 2x 10'''cm"^(Av-2). We find, however, clear deviations from this average relationship towards individual clumps. Larger than 
average A'(C'^0)/Av ratios can be found in clumps associated with the active star forming region in the northern part of Cha I. 
On the other hand, some regions in the relatively quiescent southern part of Cha I show smaller than average A'(C'**0)/Av ratios 
and also very shallow proportionality between A'(C'*0) and Ay. The shallow proportionality suggests that C'*0 is heavily 
depleted in these regions. As the degree of depletion is proportional to the gas density, these regions probably contain very 
dense, cold cores, which do not stand out in CO mappings. A comparison with the dust temperature map derived from the ISO 
data shows that the most prominent of the potentially depleted cores indeed coincides with a dust temperature minimum. It 
seems therefore feasible to use N(C^^O) and Ay data together for identifying cold, dense cores in large scale mappings. 

Key words. ISM: clouds - ISM: molecules - ISM: individual objects: Chamaeleon I - ISM: individual objects: Chamaeleon 
III 



1. Introduction 

The most common constituent of interstellar molecular clouds 
is molecular hydrogen, H2, which composes about 70 % of 
their total masses. However, direct determination of the H2 col- 
umn density is not generally possible as the molecule does not 
emit observable line radiation at low temperatures. The distri- 
bution of molecular gas, and furthermore the mass, in molecu- 
lar clouds is therefore usually derived by observing the isotopes 
of carbon monoxide, namely: '^CO, '^CO, C'^O and C"^0. 
The most common way to derive the H2 distribution is to ob- 
serve a single transition of C'^O, usually J - I -0, calculate the 
C'^O column density, A^(C'^O), by assuming the local thermo- 
dynamical equilibrium (LTE) and use a constant, "canonical", 
A^(H2) / A^(C'^O) ratio in the conversion to H2 column density. 
This method is facilitated by the facts that the transitions of 
C'^O are usually optically thin, and easily thermaUzed because 
of the relatively low dipole moment of molecule. 

There is, however, no theoretical reason for the constancy 
of the A^(H2) / A^(C"^0) ratio in the varying conditions of 
the interstellar matter In the cold and dense environment, the 
gas molecules will be depleted on the surfaces of the dust 
grains, and thus become absent from the total gaseous reservoir 

Send offprint requests to: Kainulainen, J. (jkainula@astro.helsinki.fi) 



(e.g. Little et al.T^ Burke & Hollenbach fT^S^ . This phe- 
nomenon is most obvious in dense and dark globules, in which 
abundances of molecules such as CO and CS can be decreased 
by one or two orders of magnitude towards the globule center 
(Kuiper et aL19 96 Willacy et al. ,1 998. Kramer et al. 1999, 
Alves et al. lT999l Caselli et al. lT99^ Bergin et al. 200j Jessop 
& Ward-Thompson'200P, Hotzel et al. l200Tl Hotzel et al. l200l 
Tafalla et al. 2002 1. On the other hand, molecules frozen on the 
grain surfaces can be evaporated because of heating by ener- 
getic particles or radiation. Especially, heating due to radiation 
and schocks in the vicinity of vigorous star formation should 
alter the level of depletion (e.g. Williams [19851 Charnley et 
al.il988 1. Williams ( I1985> concludes that variations can occur 
from cloud to cloud, and even within the different regions of 
the same cloud reflecting the local star formation activity. 

The dust component of interstellar matter provides another 
type of mass tracer in molecular clouds, as the dust particles 
absorb and scatter the light of the stars behind the cloud. The 
extinction of light, which is directly proportional to interven- 
ing dust column density, is very high at optical wavelengths 
and the dense cloud cores are totally opaque. However, in the 
near-infrared the extinction is only about one tenth of that in 
the visible and cloud cores are more transparent. It is therefore 
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viable to derive A^(H2) and the overall mass distribution in a 
molecular cloud with the aid of near-infrared observations. 

The A^(C'**0) / Ay ratio has been determined for several 
clouds, beginning with the work of Encrenaz et al. J1975> in 
Rho Ophiuchi. Since then, the ratio has been examined by 
Frerking et al. J1982> in Taurus and p Ophiuchi, Harjunpaa et 
al. ( 19961 in Chamaeleon I, R Corona Australis and Coalsack, 
Hayakawa et al. ( 119991 1200 lb in the Chamaeleon complex, 
and by Lada et al. <1994> . Alves et al. ( 11999 ). Juvela et al. 
( fT997l 1200 2 1 and Harjunpaa et al. (l2004t in several small glob- 
ules. These data, reviewed and compared in more detail by 
Harjunpaa et al. (I2004t . show some evidence supporting the 
correlation of A^(C'^O) / Ay with star forming activity. 

We present the derivation of the A^(C"^0) / Ay ratio in 
Chamaeleon I, and in selected region of Chamaeleon III, with 
emphasis on the possible difference between star forming and 
non-star forming regions. The present study is both deeper 
in extinction and more sensitive in C'^O than the previous 
large scale studies conducted in Chamaeleon region . In |2l 
the Chamaeleon cloud complex is briefly introduced. In ^we 
present the CO observations and NIR data used in the study. 
The relations between the derived quantities, A^(C'^O) and Ay, 
are given in |3 and these results are discussed in IS] In ^we 
summarize the results. 

2. Chamaeleon cloud complex 

The Chamaeleon cloud complex is an extensive concentration 
of molecular clouds covering dozens of square degrees on the 
sky. The complex has been studied at several wavelengths and 
it is proven very useful for studying low-mass star formation 
due to its close proximity (150 pc from the Hipparcos satellite 
data. Knude & H0g 1998.1, somewhat isolated location below 
the galactic plane {b a -16°), and the undisturbed view along 
the line of sight. Among the three large clouds in the complex, 
Chamaeleon I (Cha I) is the most active stellar factory, with 
about 200 identified Young Stellar Objects (YSOs) within its 
boundaries. In contrast, Chamaeleon III (Cha III) is a quiescent 
cloud, without current evidence for star formation. This variety 
makes the complex particularly interesting. Fig. ^ shows the 
IRAS 100 fj.m emission map of the complex, with the regions 
investigated in this study denoted. 

On a large scale, Cha I has not been mapped at many molec- 
ular line frequencies. Observations have been made in the emis- 
sion fines of H2CO and OH (T oriseva et al. 1985), '^CO, '^CO 
and C"^0 (Tori seva e t al. nM?! Mizuno et al. [T^ IT999 2001 
Haikala et al. l2005t . Tfie infrared survey using tfie DENIS 
(DEep Near-Infrared Survey of tfie soutfiern sky) pfiotome- 
try was done towards Cfia I by Cambresy et al. ([1998 ) to re- 
veal young stellar objects. Gomez & Kenyon (20011 also used 
DENIS to searcfi for YSOs and to derive tfie extinction law in 
Cfia I. Extinction maps of Cfia I fiave been previously derived 
by Toriseva et al. ( 1985) and Cambresy et al. ( 1997 1. Several 
studies fiave used different tecfiniques and instruments includ- 
ing tfie Hubble Space Telescope, tfie Very Large Telescope and 
tfie Infrared Space Observatory (ISO) to identify and cfiarac- 
terize individual objects on a smaller spatial scale (e.g. Oasa, 
Tamura and Sugitani .1999. Lefitinen et al. 120011 Persi et al. 



I2OOTI Kenyon & Gomez I2OOTI Carpenter l200l Stelzer et al. 
2004). 

Tfie previous large scale studies of Cfiamaeleon I (Mizuno 
et al.'1999 and Totfi et al. 2000 Haikala et al. 2005 1 all use 
different naming convention for tfie clumps or cores identified 
in tfie cloud. Tfirougfiout tfiis paper we use tfie naming conven- 
tion of Haikala et al. (2005) wfienever applicable. See tfie Fig. 
|2]for tfie positions of some of tfiese cores. 

Cfia III is similar in size to Cfia I but fias a lower col- 
umn density: tfie IRAS 100 jum maximum surface brigfitness 
is only about 20 % of tfiat of Cfia I. Cfia III is an example 
of a filamentary structure witfi a possible contribution of mag- 
netic fields creating a wavy, serpentine-like arcfiitecture (Gafim 
et al.[20^. Tfie cloud fias been observed in "CO and C'^O 
by Mizuno et al. (1999"), and in pait by Gafim et al. ( I2002t . 
Following tfie labelling of Gafim et al. (2 002j . an area de- 
noted as Cfialll-B fias been selected for tfiis study. It repre- 
sents one filament of Cfia III and is prominent in tfie IRAS 
A/100 = /(lOOjum) - 5/(60yum) map, containing probably tfie 
most dense matter in Cfia III. Tfie quantity A/100 has been found 
to fiave a very good correlation witfi '^CO, representing tfie 
cold dust emission in molecular clouds (Laureijs et al. 1995). 
Cfia III-B, among otfier filaments, fias been studied in more de- 
tail by e.g. Gafim et al. ( I2002t 

Tfie A^(C'^0)/Ay ratio fias been examined in tfie 
Cfiamaeleon clouds by Hayakawa et al. ( 119991 1200 It witfi 
J-band pfiotometric data from DENIS and C'**0 observa- 
tions from tfie NANTEN radiotelescope (2.7' FWHM beam). 
Hayakawa et al. (2001 1 studied tfie ratio also in different 
parts witfiin Cfia I. Harjunpaa and Mattila (1996^ derived tfie 
A^(C'**0) / Ay ratio for about 20 individual lines of sigfit across 
Cfia I using tfie SEST telescope and near-IR extintion values 
from literature. Tfie results of tfie present study are compared 
witfi tfiese results in section|5l 

3. Observations 

3.1. CO data 

Tfie C^^O (J=l-0) observations of Cfia I and Cfia III at 109 
GHz were made during several observing runs between 1988 
and 1996 witfi tfie SEST telescope on La Silla in Cfiile. Tfie 
data fiave been publisfied by Gafim et al. (2002 Cfia III), and 
by Haikala et al. ( 120051 Cfia I, H05 fiereafter). Tfie observations 
and tfie instrumentation are described in detail in tfiese articles. 
Tfie Cfia I map covers completely tfie fiigfi column density re- 
gion in tfie cloud, about 30' x 90' or 0.6 square degrees in size, 
witfi a r spacing between tfie observed positions. Tfie SEST 
beam size at tfiis frequency is 45". Tfie Cfia III map used in 
tfie present paper covers a filament called Cfia III-B in Gafim et 
al. (I2002> . witfi dimensions 10' x 30'. Tfie map spacing in tfie 
observations alternates between 30" and 1'. We fiave used tfie 
original data provided by tfie autfiors of tfie two C'**0 papers. 
Tfie data fiave been smootfied to correspond to tfie resolution 2' 
acfiieved in tfie extinction maps from tfie 2MASS data (section 

Tfie maps of integrated line area of Cfia I and Cfia III-B 
are sfiown in Fig.|2] Tfirougfiout tfiis paper, tfie C'^O emission 
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line data are presented on the radiation temperature, Tr, scale. 
To make the transformation the antenna temperatures, T^, were 
divided by the source-beam coupling efficiency 77 - 0.8, which 
is in the midway of the main beam efficiency and the moon 
efficiency of the SEST at this frequency. The intensity map 
of Cha I reveals extended features, standing out both spatially 
and in velocity. H05 used the automatic routine 'Clumpfind' of 
Williams et al. (I1994> to find out how these extended features 
decompose into numerous clumps, and studied their properties 
in more detail. The clumps with the highest peak intensities 
were found in the northern arm, close to the Cederblad 112 
reflection nebula, with masses of approximately 5 Mq. High 
intensity clumps were also found in the central region, and in 
the southern (and southwestern) parts where the C'^O emis- 
sion shows a large local maximum. Some of the most intense 
clumps are marked with numbers in the Fig.|2] 

The C'^O column densities along the sightlines were cal- 
culated from the integrated intensities using the approximation 
of optically thin line. The excitation temperature in the clouds 
was assumed to be 10 K, as was assumed by H05 and Gahm et 
al. (2002 1. The resulting C'^O column densities reach the max- 
imum values 2 x 10'^ cm"^ and 2 x 10'^ cm"^ in Cha I and Cha 
III-B, respectively. 

3.2. NIR data 

We derive the extinction maps of Cha I and Cha III-B using 
JHK near-infrared photometry of 2MASS (2 Micron All-Sky 
Survey). The 10-crlimitting magnitudes of the 2MASS data are 
15.8, 15.1 and 14.2 in J, H and K bands, respectively. In par- 
ticular, we use the NICER color excess method (Near-Infrared 
Color Excess Revised) presented by Lombardi & Alves (2001 1, 
which generalizes the color excess method formulated in Lada 
et al. ([1994 1 for multi-band observations (see Lombardi & 
Alves|2001 for the complete description of the method). The 
data of the 2MASS survey in three bands yield two indepen- 
dent colors, namely J — H and H - K, which are used in the 
following to implement the NICER method. Before any other 
calculations, the 2MASS photometry is tranformed to the CIT 
(Caltech Institute of Technology) photometry using the trans- 
formation factors given in Eqs. (12)-(15) of Carpenter (2001 1. 
The color excesses of background stars are scaled to corre- 
spond to the extinction in the V band using the diffuse dust 
reddening law (Ry = 3.1) of Mathis ( fT990l : 

Av = 15.87 X £H-/r (1) 

Av = 9.35 X Ej^H (2) 

It should be noted that the calculated "visual extinction" actu- 
ally is a scaled version of the infrared extinction and does not 
reflect variations due to possibly altering value of Ry- 

The 2MASS photometry of two on-source frames and two 
reference frames was obtained from the 2MASS All-Sky Point 
Source Catalog ' . The two on-source fields covered the cloud 
areas of interest, Cha I and Cha III-B, and they contained ap- 
proximately 13000 and 3000 stars, respectively. The complete 

' http://irsa.ipac.caltech.edu/applications/Gator/ 
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Table 1. New sources with NIR excess emission in Cha I 



Star 


a (2000.0) 


S (2000.0) 


J-H 


H-K 


1 


llh04m 50.8s 


-76d51m20s 


0.77 


0.74 


2 


llh04m 59.6s 


-77d 35m 58s 


0.82 


0.81 


3 


llh07m 19.3s 


-76d 54m 19s 


0.63 


0.65 


4 


llh 10m 57.5s 


-77d 44m 44s 


0.75 


0.69 



lists of sources in these fields are referred to as the full samples 
in the following. Two reference fields were needed to calculate 
the intrinsic colors of stars that are free from extinction. The 
reference fields were selected using the IRAS 100 fj.m emis- 
sion map of the complex as a guideline, and circular areas with 
radii of 10 arcminutes coiTesponding to the minima of 100 /urn 
emission were chosen. The locations of these four fields are 
outlined in Fig.^ 

3.2.1 . Contribution of young stars 

About 200 young stellar objects have been identified within 
the boundaries of Cha I in previous studies. Due to the sur- 
rounding dust envelope these stars generally show NIR-excess 
emission with respect to the main sequence stars, and thus the 
extinction value derived towards them is incorrect. Especially 
in the densest regions, which are likely to contain YSOs, the 
extinction estimation relies on the photometry of a relatively 
low number of stars. If one (or several) of these stars possess 
abnormal colors, the calculated extinction may differ greatly 
from the "real" value. It is therefore important, that the young 
stars are removed from the sample, even though their number 
is small compared to the full sample. 

To identify the young stars in Cha I we used the list of 
known YSOs from Carpenter GOOD complemented with the 
sources from Gomez & Kenyon (2001). The 2MASS sources 
with positions and colors matching with known YSOs were re- 
moved from the sample. Altogether ~ 90 % of the young stars 
were recovered and removed from the full sample (positional 
matching). The remaining set of stars, from which the YSOs 
have been removed, is referred to as the restricted sample in 
the following. 

The color-color diagrams of 2MASS sources having pho- 
tometry with S /N > 10 are presented in Fig.|3 with intrinsic 
color-curves of unreddened stars superposed (Bessell & Brett 
1988). The reddening sector is indicated with dashed lines, and 
it corresponds to the reddening law Ej^u/Eh-k = 1.8 + 0.03 
derived in Cha I by Gomez & Kenyon ( 1200 U . The previously 
known YSOs are marked in the diagram with diamonds, and 
obviously no new candidates with high excess emission can be 
claimed. A few sources fall to the right from the reddening sec- 
tor around {H - K,J - H) k (0.8, 0.7). It is possible that these 
sources are young stellar objects, even though the colors imply 
small extinctions and thus small column densities of surround- 
ing dust. Their locations are given in Tab.^ 
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3.2.2. Visual extinctions 

Using the restricted sample, from which the contribution of 
young stars has been eliminated, the NICER color excess 
method is applied. The extinction in the densest regions of Cha 
I rises high enough to result in regions of a few arcminutes in 
size where no 2MASS sources are detected. By requiring at 
least one star within the FWHM range from the center of each 
pixel, and estimating the surface density of stars from the refer- 
ence field, we expect that the extinction estimation is limited to 
Ay < 20'" on the average. This, of course, varies according to 
what is the true stellar surface density along each sightline. To 
achieve reasonable spatial resolution over most of the Cha I we 
select the FWHM resolution of the NICER method to be 2.0', 
even though it slightly falls below the diameter of some star- 
less regions. The requirement of having at least one star within 
the FWHM range is thus not met for some pixels, and the de- 
termination of extinction in those pixels is impossible. In Cha 
III-B where the column density is generally smaller we select 
the FWHM resolution to be 1.5'. The resulting extinction maps 
are presented in Figs.0]and|5l 

The errors of the derived extinction maps originate from the 
uncertainty of the intrinsic colors, (J - H)q and (H - K)o, and 
the uncertainty of the observed magnitudes of the on-field stars 
(Eqs. 9 and 16 in Lombardi & Alves l200Tl . Generally, the 1-cr 
relative error over the maps is about 6%. On average, the error 
of extinction per pixel is expected to increase with decreasing 
star density in denser regions, as the photometric errors become 
larger. Fig. |6l describes the behaviour of the average standard 
deviation of pixels as a function of the number of stars within 
the FWHM radius. In the low column density regions there are 
over twenty stars within the FWHM range and the average stan- 
dard deviation is about 0.2'". When the star density drops, the 
average standard deviation increases and reaches about 0.5'" at 
point when there is only one star within the FWHM range. The 
maximum eiTor in this case is about 1.5'", which equals to the 
relative error of about 8 %. Thus, even in the case where the 
extinction measurement relies on a small number of stars, we 
find the error acceptable. In Cha III-B the extinction is gen- 
erally lower and the probability of having apparently brighter 
stars within the FWHM range is larger. This means that the av- 
erage photometric errors are lower than in the region of higher 
extinction. Thus, the increase of the average error at low star 
densities is not as steep as in Cha I. 

The extinction maximum in Cha I is located in the 
Cederblad 110 region (central part), and it rises up to 19'" 
at the selected resolution. As mentioned above, extinction in 
some pixels remains undefined as no stars lie within the ra- 
dius defined by the FWHM resolution (these "empty" pixels 
are marked white in the extinction map). The extinction in these 
regions is higher than can be determined. In the Cederblad 112 
region (northern part) the maximum extinction is about 15'", 
and the "arm" southwards from Cederblad 112 has Ay > 10'". 
Numerous clumps with Ay > 10'" can be found all over the 
cloud, even in the more quiescent southern part. A large num- 
ber of less opaque clumps with Ay a 5'" are present on the 
outskirts of the cloud, being most clearly visible west of the 
northern arm. 



All four C'^O cores detected by Mizuno et al. J1999> lo- 
cated within the borders of the extinction map are clearly seen 
as local extinction maxima. The southernmost of these appears 
to have very high column density as the 2MASS source den- 
sity decreases to zero. In addition, most of the very cold cores 
(VCCs) detected in the ISOPHOT Serendipity Survey by Toth 
et al. (2000) in this region are recovered as well. See Fig.l^for 
the positions of these cores. 

In H05 the central and southern regions of Cha I are re- 
vealed as concentrations of several clumps, some of which re- 
main unresolved at the resolution of the extinction map. Even 
though the clumps are not individually resolved, their inte- 
grated effect is visible. The most obvious concentrations de- 
tected are clumps 1 and 2 in the north, 4 and 7 in the east, 5 and 
6 in the south and 8 (among others) in the center. Also many 
less intense clumps, such as 10, 12 and 21, are particularly well 
resolved. As a bottom line, even most of the low intensity (and 
less massive) C'^O clumps such as 35, 38 and 53 seem to have 
a well defined counterparts in the Ay map. 

Cha III-B has the maximum extinction of about 10"' occur- 
ing in the southern part of the studied area. The region of the 
highest extinction breaks into two maxima, which is also seen 
(more clearly due to better resolution) in the Fig. 7 of Gahm et 
al. ('2002V Many smaller features of the C'**0 map such as the 
local maximum in the north and the short arm leading south- 
east from the southern maximum are also visible, but somewhat 
smoothed by the lower resolution of the extinction map. 

4. Correlations between N(C^^O) and Ay 

After smoothing the Ay and A^(C'**0) maps of Cha I and Cha 
III-B to the common resolution of 2.0', the C"*0 column den- 
sities and visual extinctions are compared on a pixel-by-pixel 
basis. The resulting coiTelation diagrams for Cha I and Cha III- 
B are are shown in Figs.^t ™d b. A linear fit to all data points 
in Cha I gives the following relationship: 

N(C^^O) = (1.9 + 0.1) X lO^^Ay 

- (4.9 + 0.7) X 10'"^ (Chal) (3) 

The plot in Fig.^ suggests that the proportionality is steeper 
than the fit at low values of Ay, but becomes slightly shallower 
above Ay > 7'", and finally at Ay > 12'" most datapoints lie 
below the fitted line. 

In Cha III-B we obtain the following fit: 

A?(C'^0) = (2.0 ± 0.1) X 10'**Ay 

- (3.9 + 0.4) X lO''* (ChalllB) (4) 

In Cha III-B the extinction remains below 10'", and no clear 
deviation from the linear relationship can be seen. The parame- 
ters of the linear fits to the Cha I and Cha III-B data sets do not 
differ significantly at this scale, where the relation is averaged 
over the Cha I cloud. 

The extinction map of Cha I has several opaque spots 
with Ay > 10'". Almost all of them can be identified with 
C'^O clumps found by H05. These clumps with radii ranging 
from about 4' to 10' can be discerned also on the C'**0 map 
smoothed to 2'. The A'(C'^O) vs. Ay plots of the regions of 
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seven selected clumps are shown in Figs.Q;-i. Figs0; and i rep- 
resent regions where the A^(C'^0)/Av ratio is higher, whereas 
Figs-Et ™d h represent regions where the A^(C'**0)/Av ratio is 
lower than expected from the Unear fit to the whole Cha I. 

5. Discussion 

As the CO molecule is chemically very stable, the variations 
of the A^(C"*0)/Av ratio depend on the accretion of molecules 
onto grain surfaces and on desorption processes which release 
them back into the gas phase. In the obscured and dense cores 
with no central source the abundance drop of CO, caused by 
accretion onto the surfaces of cold dust grains, may be several 
orders of magnitudes. The dust temperature is very low, about 
~ 10 K and evaporation due to thermal desorption very unefli- 
cient. The rate of the thermal desorption can be written as: 

^(M)c.exp{-^} (5) 

where Eb(M) is the binding energy, the dust temperature 
and M refers to the composition of the grain mantle (Watson 
& Salpeter [l972> . Non-thermal processes such as desorption 
caused by cosmic rays are also present, and in fact required 
to sustain the observed abundances in cold cores (Willacy and 
Williams 1993 Hasagawa and Herbst 1993 1. However, in the 
presence of a central heating source the thermal desorption of 
CO quickly overcomes other processes due to the low binding 
energy of the CO molecule. In the envelope surrounding the 
heating source the CO is evaporated in a short timescale, and 
most of the CO should remain in gas phase as the dust temper- 
ature keeps rising (e.g. Rodgers & Charnlev 120031 Lee et al. 

The depletion of gas phase C'^0 can be quantified with the 
aid of the depletion factor, /d, which is defined as the ratio 
of the total fractional abundance to the gas phase abundance: 
/d = Af(C'**0, gas + dust)/;(f(C"^0, gas). As discussed in Hotzel 
et al. J2002I I. in an isothermal cloud where these processes are 
in equilibrium /d is proportional to the gas density according 
to /d = 1 + CnHj, where C = A/B is the ratio of the adsorption 
and desorption coefficients. These coefficients depend on the 
grain properties, temperature and the flux of heating particles 
(primarily cosmic rays). A reasonable set of parameters gives 
the following estimate in dark clouds: C ~ 5 10 "^cm^ (Hotzel 
et al. 1200 2). Large changes in adsorption and desorption coef- 
ficients are not likely to occur within a molecular cloud, unless 
local heating by radiation or shocks increases the grain temper- 
ature above a critical temperature of CO desorption, which lies 
in the range 20-30 K (Collings et al. 2003). The no-depletion 
fractional C"^0 abundance, ;^'(C'^ (9, gas + dust), is determined 
by the gas-phase production of CO in the low-density cloud en- 
velope, and is probably constant throughout the cloud. 

The density dependence of the depletion factor is trans- 
ferred also to the observed A^(C'^0)/Av ratio: 

A?(Ci^O,gas) _ J n(C^^O, gas) ds 

Av a Jn(ll2)ds 

J n(Ci^O,gas + dust)/(l + Cn(}i2))ds 

a J nQi2)ds 



/^(C'«0,gas + dust)^^^. 

- p ■ (o) 

a J n(H2)c/i 

The Av to A^(H2) conversion factor, a, is given in Eq. ^2 For 
example, assuming line-of-sight homogeneity one obtains the 
following relation: 

N(C^^O, gas) _ ;t'(C"^0, gas + dust) 1 

M ~ a 1 +Cn(H2) ■ 

Although this assumption is not very realistic, the example il- 
lustrates the tendency that the measured A^(C"*0)/Av ratios are 
likely to be smaller along line of sights which contain very 
dense gas than elsewhere. 

Based on the arguments above, the variations of the 
A^(C'^O) / Av ratio can be attributed to two main processes 
which can be seen as consequencies of gravitational contrac- 
tion. On the one hand, a large column density required for C'^O 
detection shelters C'^O efficiently from the external radiation 
field, and as the accretion rate of molecules onto the grain sur- 
faces depends strongly on the density, a high degree of CO de- 
pletion can be expected towards dense cores without internal 
heating sources. On the other hand, the densest cores are likely 
to be active sites of star formation which can give rise to an 
increased local heating releasing molecules from grains. Also 
violent events such as winds and jets from newly born stars 
may play a role in this process. 

Apart from the physical reasons, the validity of the assump- 
tions made in the calculation of CO column densities may vary 
as a function of density (i.e. optical depth). The comparison 
of C'**0 and C'^O line profiles in H05, however, supports the 
view that the C'^O line remains optically thin. The assump- 
tion of the constant excitation temperature (10 K) is perhaps 
more fragile in this sense, as the temperatures in cloud cores 
have been predicted to drop towards the core center (Zucconi 
et al. 12001; . Galli et al. I2002> and observations have also given 
such indications (e.g. Ward-Thompson et al. 2002 Lehtinen et 
al. 2004 ). The observed temperature gradients are nevertheless 
quite small, and even going from, say, Tgx = 8 K to Tgx =15 
K does not produce a large enough difference to explain the 
observed N(C'*'0) variations in clumps (see the discussion e.g. 
in Harjunpaa and Mattila 1996). Furthermore, the neglection 
of the temperature gradient should lead to similar variations in 
clumps with equal maximum extinction. The situation in Figs. 
Eh-i representing clumps with equal maximum extinctions sug- 
gests therefore that temperature gradients cannot be the main 
reason for the variation in the A^(C'^O) / Av ratio. 

5.1. N(C^^O) / Av ratio in Cha I and III-B 

The A^(C"*0) /Ay relations derived for Cha I and Cha III-B do 
not differ significantly from each other on the large scale where 
the whole Chamaeleon I region is considered (Eqs.|3jE}. This 
is not in agreement with the result derived in Hayakawa et al. 
J200H . in which a significant difference between the clouds 
was found. In Hayakawa et al. ( 2001 ) the slope derived for Cha 
I was steeper than that of Cha III by a factor of 1 .4. The rela- 
tionships derived by Hayakawa et al. (1200 1> are quoted below: 

A^(C'**0) = (3.1 +0.1)x lO'^^Av 



6 



Kainulainen, J. et al.: The ratio of A'(C'*0) and Ay in Chamaeleon I and III-B 



- (4.4 + 0.7) X 10' 



(Chal.Hayakawa) (8) 5.3. Variations ofN(C^^O) /Ay in Clia I 



NiC^^O) = (2.3 + 0.3) X lO'^^Av 
-(0.3 + 0.8)x lO"* 



(Chain, Hayakawa) (9) 



A possible reason for this difference in results appears to 
be discrepancy between the C'^O intensities observed with the 
SEST and the NANTEN telescopes: the intensities observed 
using SEST in Cha I are smaller by a factor of 1.5-2.2 com- 
pared to intensities at the positions of cores given in Mizuno et 
al. 



Ill999> . On the other hand, the intensities of Mizuno et al. 
<1999> in Cha 111-B exceed those of the present study only by a 
factor of ~ 1.2. The comparison was done after smoothing the 
SEST maps to correspond the 2.7 FWHM beam of NANTEN. 
However, the number of common datapoints especially in Cha 
III-B is small, and a proper overall comparison is not possible. 
But if the difference in the intensities is generally at this level, it 
is alone sufficient to explain the difference between the results 
of the present study and Hayakawa et al. (.2001 J . 



5.2. NiH2) /N(C^^O) ratio and tine cloud masses 

To express the observed A^(C'^O) / Ay relation in terms of hy- 
drogen column density we use the A^(HI-hH2) / ^(B-V) ratio 
derived from the observations of the Copernicus satellite by 
Bohlin et al. dflWSb: 



N(H + H2) 



= 5.8 X 10^'cm^^mag 



(10) 



E{B - V) 

The sightlines of Bohlin et al. ( 11978 1) were restricted to low 
extinctions (Ay < 2), but the result has been verified also for 
denser clouds (Vuong et al. 2003). Adopting the value Ry = 
Ay /Es-v = 3.1, consistent with the implementation in NICER, 
we have: 

^(H2) 



= 9.4 X 10^"cm"^mag 



(11) 



We first calculate the masses of the clouds from the extinc- 
tion maps using Eq. The masses are 260 Mq and 18 Mq 
for Cha I and Cha III-B, respectively. Then, we apply Ea.llOlto 
Eqs.|3land|3 and have: 



A^(H2) = 4.9 X 10*'A^(C''*O) + 2.4 x 10^'cm" 



(Chal) (12) 



A?(H2) =4.7xl0''A^(C"*O)-Hl.8xl02'cm"2 (ChaIIIB)(13) 
These relations are actually very close to the "canonical" C'^O 



/ Ay values derived by Frerking et al. (I1982I I: 

A^(H2) = 5.9 X lO^NiC^^O) + 1.3 x lO^'cm"^ (Taurus)(14) 

A^(H2) = 5.9xl0''A^(C"*O)+3.9xl0^'cm-2 (pOphiuchi)(15) 

The masses of the observed areas using Eas.ll2landll3lare 
280 Mq and 21 Mq. The masses derived using the two different 
approaches are not independent and turn out to be, as expected, 
close to each other. This reflects the fact that most of the mass 
comes from the regions with moderate density, that is, from the 
region where the linear A^(C'^O) / Ay relation holds well on the 
average. 



Figs-Et-i point out differences in the A^(C'**0) vs. Ay relation 
when comparing individual clumps in Cha I. The ratio in the 
northern region differs significantly from the rest of the cloud, 
and small differences between other clumps can be seen mainly 
in the overall C'^O level (the intercept of the linear model). 
Also, few interesting exceptions are found, especially in the 
western and southern parts of the cloud. 

The combination of clumps 1 and 2 in the north represents 
the high-CO side of the A^(C"*0) vs. Ay relation over the full 
Ay scale (Fig.Et)- Compared with clumps in less active regions 
the average A^(C'^O) level impUes the overabundance of C'^O 
by a factor of 1.3 + 0.2 with respect to the relation averaged 
over the whole Cha I, and by a factor of 3.5 ± 0.5 with respect 
to the relation of the clump with low C'**0 level in Fig. Eh- 
This supports the view that in this region the C'^O abundance is 
indeed enhanced due to very high star formation activity. Along 
the "arm" southwards from clumps 1 and 2 there is another very 
massive clump, clump 12, which has no identified YSOs in it. 
The extinction rises close to the value of the clumps 1 and 2, 

1. e. ~ 12"', but the C'^O level remains significantly lower. The 
plot shown in Fig. Efe is different from that of clumps 1 and 

2. Again, this tendency could be interpreted in terms of less 
intense radiation, and thus higher overall level of gas depletion. 

The active region in the central part of Cha I has a very high 
column density with Ay > 15"' over a large area. For several 
tens of pixels the star density drops to zero per NICER beam, 
i.e. the FWHM range of NICER, and the extinction cannot be 
determined. Also, relying on very few stars per beam, the error 
in the pixels where the measurement is possible becomes large 
(see Fig.|6|l. The datapoints spread out over a wide range in the 
A^(C'^O) vs. Ay correlation diagram (Fig. 01)- The linear fit to 
the datapoints follows the Eq.|3very closely. 

The westernmost C'^O maximum including clumps 4 and 
7 shows a peculiar behaviour with close to constant A^(C'**0) 
level over the range Ay ^ 5-14"'. This results in a narrow, 
almost horizontal relation implying that the level of depletion 
in this clump is particularly high (Fig. Et)- A few identified 
YSOs lie loosely scattered in this region, but generally it be- 
longs to the low-activity part of the complex. The high level 
of depletion, however, is an implication of large density and 
low temperature. Fig. [S] showing a tentative dust temperature 
map derived from the ISO 100, 150 and 200 jim observations 
(Lehtinen et al. .2005.. in preparation) indeed points out that this 
region represents the coldest region of the southern part of the 
cloud. We thus identify the region as a prominent, cold, heavily 
depleted dense core. 

The southernmost of the large C'^O maxima, consisting 
of clumps 5 and 6 suffers from the same lack of stars due to 
high extinction as the central region. The maximum extinction 
derived is about 14"', but taking the region without any back- 
ground stars into account, the maximum extinction is substan- 
tially higher making this region most likely the densest part of 
the southern side of the Cha I cloud. This clump also corre- 
sponds to the position of a very cold core detected by Toth et 
al. (.2000j from the ISOPHOT Serendipity Survey data, having 
an ISO /(170)//(100) color temperature < 12.5 K. The Fig.|8] 
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confirms that the clumps indeed coincide with the local tem- 
perature minimum. 

For comparison, a pair of clumps with very different 
vs. A\j relations are shown in Figs.l7h-i. Spatially the 
clumps belong to the southern part of Cha I and are separated 
only about 15'. The former one has only very weak counter- 
part in C'^O emission, even though its extinction rises above 
Ay=\(y\ The latter one, located close to a group of identified 
YSOs has a much stronger C'^O emission compared to the ex- 
tinction. 

On the basis of the earlier discussion, it seems likely that 
the clumps where the A^(C'^O) vs. Ay proportionality is shal- 
low (Figs.Et ™d h) are denser than those which roughly follow 
the average relationship (e.g. Fig.^). The latter can represent 
column density maxima without well-defined centres of grav- 
ity. On the other hand, in the two clumps in the northern part 
of Cha I (Fig.|7j;) the A^(C'**0)/Av ratio is larger that the av- 
erage. This region has recently formed a dense cluster of stars, 
and is associated with an extented bipolar outflow (Mattila et 
al. I1989I I. It seems therefore possible that the relatively large 
abundance of gaseous C"^0 in these clumps is a result of inten- 
sified desorption related to heating by shocks or radiation from 
newly born stars. 

6. Conclusions 

We have derived the maps of visual extinction in the ac- 
tive stellar star forming cloud Chamaeleon I and in a dense 
quiescent filament of the Chamaeleon III cloud using the 
NICER color excess technique and 2MASS data. We have stud- 
ied the A^(C'^O) / Ay ratio and searched for differences be- 
tween the two clouds, and between the regions and clumps in 
Chamaeleon I. For this study we have used the C'^O data from 
the SEST telescope published in Haikala et al. (2005 Cha I) 
and Gahm et al. J2002i Cha III), by smoothing them to the res- 
olution of the extinction maps. We have also used 2MASS data 
to derive the reddening law in Chamaeleon I. The conclusions 
of our work are summarized in the following: 

1 . The derived visual extinction maps of Cha I and Cha III-B 
resemble well the complex column density structure seen 
in the C'^O maps. With the 2' resolution the maximum ex- 
tinction reaches to Ay ^ 19™ in Cha I and Ay > 10'" in Cha 
III. Several opaque spots with Ay > 10"' can be found on 
the Cha I map, and most of them can be attributed to the 
C"^0 clumps identified by Haikala et al. (p005l. 

2. The average relationship between A^(C'^O) and Ay is simi- 
lar in Cha I and Cha III, i.e. N(C^^O) ~ 2 x lO^^* cm-2(Ay - 
2). This is in contrast with the previous result of Hayakawa 
et al. ( 2001 ), who find a steeper proportionality in Cha I us- 
ing NANTEN observations. The reason for this discrepancy 
is probably related to an inexplicable difference between 
the intensity scales of the two C'^0 data sets towards Cha 
I. 

3. The color-color diagram (H-K,J-H) of 2MASS stars in 
Chamaeleon I reveals four new YSO candidates with NIR 
excess emission. The new sources with excess emission are 
listed in Tabled 



4. Regional variations of the A^(C'^O) / Ay ratio are found 
within Chamaeleon I. The most obvious example is 
the northernmost active star formation region where the 
A^(C"^0) level is enhanced by a factor of 1.3 + 0.2 with 
respect to the average relation in the cloud, and by a fac- 
tor of 1.5-3.5 with respect to the clumps with C'^O level 
lower than average. We suggest that this is caused by an 
increased amount of gas phase C'^O as a result of active 
star formation. On the other hand, regions with a markedly 
shallow A^(C'^O) vs. Ay proportionality can be identified. 
These are likely to represent cold, dense cores with a high 
degree of molecular depletion. It looks thus conceivable to 
use the combinaton of 2MASS data and large scale C'^O 
surveys for finding low-temperature molecular cores. This 
suggestion is substantiated by the fact that the most exten- 
sive of the cores standing out on the A^(C"^0) vs. Ay plot 
coincides with a local dust temperature minimum discov- 
ered using ISO data. 
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Fig. 1. IRAS 100 yum emission map of the Chamaeleon cloud complex. Chamaeleon I and III-B are marked with rectangles, and 
positions of the reference fields used in the extinction calculation are marked with circles. 
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Fig. 2. Left: Integrated C'**0 intensity map of Chamaeleon I. The map is smoothed to correspond the 2.0' FWHM resolution 
of the extinction map (Fig.|3}. The contour levels start from 0.5 K km s"' and the step is 0.2 K km s"'. The numbers give the 
positions of some CO clumps identified by H05. Right: The same for Cha III-B. The contour levels start from 0.4 K km s ' and 
the step is 0. 1 K km s ' . The letters mark the most intense clumps identified by Gahm et al. J2002> . 
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Fig. 3. Left: Color-Color diagram of 2MASS stars with S/N > 10 in Cha I. Diamonds mark the previously identified YSOs from 
the list of Carpenter (2002). The solid curves are the intrinsic colors of unreddened stars (after BesseU & Brett 1988), and dashed 
lines tangenting these curves define the reddening sector. Right: The same for Cha III-B. Plus signs in the upper left comer of 
diagrams are the typical errors of colors. 
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Fig. 4. Visual extinction map of Chamaeleon I derived using the NICER color excess method (Lombardi & Alves 2001) and JHK 
photometry of 2MASS survey. The FWHM resolution of the map is 2.0'. The white pixels result from not having a single star 
within the FWHM range from the pixel centre, which means that the extinction in these pixels is higher than can be determined. 
The contour levels start from Ay = 5™ and the step is 2™. The letters mark the position of CO cores investigated in Fig. 
White circles indicate the C'^O cores of Mizuno et al. ( 1999) and the plus signs the very cold cores detected in the ISOPHOT 
Serendipity survey by Toth et al. (I2000> . 
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Fig. 5. Visual extinction map of Chamaeleon III-B derived using the NICER color excess method (Lombardi & Alves 2001) and 
JHK photometry of 2MASS survey. The FWHM resolution of the map is 1.5'. The contour levels start from Ay - 4"' and the 
step is l". 
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Fig. 6. Left: The average standard deviation of extinction as a function of star density in Cha I. The x-axis gives the number of 
stars within the FWHM range from the pixel center. Right: The same for Cha III-B. 
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Fig. 7. The correlations of A^(C'*'0) and Ay in Cha I (a), Cha III-B (b), and in selected clumps of Cha I. The clump locations 
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Fig. 8. A tentative dust temperature map of Cha I derived from the ISO 100, 150 and 200 fxm observations (Lehtinen et al. in 
prep.). The numbers refer to some of the C'^O clumps identified by H05. The temperature minimum coincides with the region 
of the clumps 4 and 7, revealing a prominent cold, heavily depleted core. 



